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Abstract. We study new types of geometric query problems defined as follows: given a geometric set 
P, preprocess it such that given a query point q, the location of the largest circle that does not contain 
any member of P, but contains q can be reported efficiently. The geometric sets wo consider for P are 
boundaries of convex and simple polygons, and point sets. While we primarily focus on circles as the 
desired shape, we also briefly discuss empty rectangles in the context of point sets. 



1 Introduction 

Largest empty space recognition is a classical problem in computational geometry, and has applications in 
several disciplines like data-mining, database management, VLSI design, to name a few. Here the problem 
is to identify an empty space of desired shape and maximum size in a given region containing a given set 
of obstacles. Given a set P of points in M^, an empty circle, is a circle that does not contain any member 
of P. An empty circle is said to be a maximal empty circle (MEG) if it is not fully contained in any other 
empty circle. Among the MECs', the one having maximum radius is the largest empty circle. The largest 
empty circle among a point set P can easily be located by using the Voronoi diagram of P in 0(n log n) time 
[25]. The maximal empty axis-parallel rectangle (MER) can be defined in a similar manner. The literature 
on recognizing the largest empty axis-parallel rectangle among obstacles has spanned over three decades in 
computational geometry. The pioneering work on this topic is by Namaad et al. [20] where it is shown that 
the number of MERs' (m) among a set of n points may be O(n^) in the worst case. In the same paper, an 
algorithm for identifying the largest MER was proposed. The worst case time complexity of that algorithm is 
0(mm(n^, mlogn)). The best known result on this problem runs in 0(nlog^ n) time in the worst case. The 
same time complexity result holds for the recognition of the largest MER among a set of arbitrary polygonal 
obstacles [21]. However, the largest MER inside an n-sided simple polygon can be identified in O(nlogn) 
time [5] . The worst case time complexity for recognizing the largest empty rectangle of arbitrary orientation 
among a set of n points is O(n^) [9]. 

Although a lot of study has been made on the empty space recognition problem, surprisingly, the query 

version of the problem has not received much attention to the best of our knowledge. The problem of finding 
the largest empty circle centered on a given query line segment has been considered in [4] . The preprocessing 
time, space and query time of the proposed algorithm are O(n^logn), O(n^) and O(logn), respectively. In 

practical applications, one may need to locate the largest empty space of a given shape in a desired location. 
For example, in the VLSI physical design, one may need to place a large circuit component in the vicinity of 

* This version is a significant update of our earlier arXiv submission arXiv:1004.0558vl. Apart from new variants 
studied in Sections 3 and 4, the results have been improved in Section 5. Part of this work was done when the 
authors were visiting IISc, Bangalore, and TIFR, Mumbai. 



some already plaeed components. Such problems arise in mining large data sets as well, where the objective 
is to quickly study the characteristics (such as the area of the empty space) near a query point. 

In this paper, we will study the query versions of the empty space recognition problem. If the desired object 
is a circle, the problem is referred to as maximal empty circle query (QMEC) problem, and if the desired 
object is an axis-parallel rectangle, the problem is referred to as maximal empty rectangle query (QMER) 
problem. The following variations are considered. 

Given a convex polygon P, preprocess it such that given a query point q, the largest circle inside P that 

contains the query point q can be identified efficiently. 
Given a simple polygon P, preprocess it such that given a query point q, the largest circle inside P that 

contains the query point q can be identified efficiently. 
Given a set of points P, preprocess it such that given a query point g, the largest circle that does not 

contain any member of P, but contains the query point q can be identified efficiently. 
Given a set of points P, preprocess it such that given a query point q, the largest rectangle that does not 

contain any member of P, but contains the query point q can be identified efficiently. 

We believe that our work motivates study of new types of geometric query problems and may lead to a very 
active research area. The main theme of our work is to mainly understand which problems can be solved 
in subcjuadratic preprocessing time and space, while ensuring polylogarithmic query times. Our results are 
summarized in Table 1. 

Table 1. Complexity results of different variations of largest empty space query problem 



Geometric set 


Shape of 
empty space 


Preprocessing 
time 


Space 


Query time 


Sections 


Convex Polygon 


circle 


0{n) 


0{n) 


O(logn) 


3 


Simple Polygon 


circle 


0{n log'^ n) 


0{n log'^ n) 


0(log^T7,) 


4 


Point Set 


cirek^ 


0[,r\og,t) 




O(log;0 


■5 


Point Set 


rectangle 


0{n'^ logn) 


0{ri^ log n) 


O(logn) 
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In the course of studying these problems, we developed two different ways of implementing a key data 
structures for storing n circles of arbitrary sizes such that when a query point q is given, it can report the 
largest of the circles that contains q. This data structure may be of independent interest since it may aid in 
several other geometric search problems. 

2 Preliminaries: LCQ-problem 

In this section, we want to build a data structure called the largest circle query data structure (LCQ. in short) 
for the point location in an arrangement of circles. Our input is a set C = {Ci, C2, . . . , C„} of circles in 
non-increasing order of their radii. In the preprocessing phase we will construct the data structure. When a 

query point q is given, it must report the largest circle in C that contains q 01 a. null value if q is not enclosed 
by any circle in C. For simplicity, we assume that at most two circles intersect at any point on the plane. 

We provide two ways of building the LCQ data structure. The first method uses divide-and-conquer, leading 
to a solution that is optimized for preprocessing time. The second method uses a line sweeping technique 
similar to [24] , and it gives a solution with better query time. The complexity results are given in Table 2 

2.1 A divide-and-conquer solution 

Preprocessing: We form a tree V of depth O(logn) as follows. Its root r represents all the members in 
C, and is attached with a data structure vor(r) with the circles in C(r) = C. The two children of root, say 
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Table 2. Table comparing the two solutions for LCQ 



Techniques 


Preprocessing time 


Space 


Query time 


Divide- and-conquer 


0(n log^ n) 


O(nlogn) 


O(log'n) 


Line sweep 


0{n^ logn) 


0{n') 


0(log n) 



Vi and 2? J., represent the sets C(r^) = {Ci, C2, . . . , C'|^ } and C{r^) = {C|^|j+i, . . . , C„_i, C„}, respectively. 
These define the associated structures vor(r£) and vor(rr-) of and Vj. respectively. The subtrees of Vg, and 
T>r are defined recursively in the similar manner. Finally, the leaves of T> contain Ci, C2, . . . , C„, respectively. 
The tree is computed in a bottom-up fashion starting from the leaves. The task of the data structure vor(u) 
associated to a node v is to efiiciently report whether or not the query point q lies inside the union of the 
circles in C{y) it represents. We will use Voronoi diagram in Laguerre geometry of the circles in C{v) [15]. 
Each cell of this Voronoi diagram is a convex polygon and is associated with a circle in C{v). The membership 
query is answered by performing a point location in the associated planar subdivision. For a node v, vor(u) 
can be computed in 0(|C(w)| log |C(w)|) time and membership query can be answered in 0(log |C(w)|) time 
[15]. 

Query answering: To find the largest circle in Cq containing the given query point q, we start searching 
from the root r of V. If q does not lie in the union of circles C(r) = C, then q is contained in an empty 
circle of size infinity. We need not proceed further in the tree. However, if the search succeeds, we need to 
continue the search among its children. A successful search at a node 77 indicates that q must lie cither in 
the union of circles of its left child or the right child or both. We first consider its left child ve, that contains 
the larger l^Mi circles of node v. We search in the associated structure voT{ve). If the search succeeds (i.e., 
q e ^cecivf:)^^)^ the search proceeds in the subtree rooted at V£. However, if the search fails, surely q lies in 
the union of circles C{vr), and the search proceeds in the subtree rooted at Vr- Proceeding similarly, one can 
identify the largest circle Cq containing the query point q. 

Theorem 1. A set C of n circles can be preprocessed in O(nlog^n) time and 0{n\ogn) space so that LCQ 
queries can be answered in O(log^n) time. 

2.2 A line sweep solution 

We assume a pair of orthogonal lines on the plane to represent the coordinate system. The circles arc given 
as a set of tuples; each tuple representing a circle consists of the coordinates of its center and the radius 
of that circle. The ordered set of vertices V = {vi,V2, ■ ■ ■ ,v\v\) of the arrangement A{C) consists of the 
(i) leftmost and rightmost point of each circle in C, and (ii) points in which a pair of circles intersect. We 
assume, further, that the vertices of A{C) have unique x coordinates, thereby allowing the elements of V 
to be stored in increasing order of their x coordinates. A maximal segment of any circle in C that does not 
contain a vertex is called an edge of A{C). We use E to denote the set of all edges of A{C). Since we include 
the left and right extremeties of a circle in the set of vertices, the edges are always x-monotone. Each edge 
is attached with two fields IDi and ID2 indicating the largest circle containing the cell above and below it 
respectively. Each cell is bounded by the edges of A{C), and is attached with an index ID indicating the 
largest circle containing that cell. We compute the arrangement A{C) as follows: 

Step-1 Cut each circles into pseudo-segments such that a pair of pseudo-segments intersect in at most one 
point. If any of these segment contains the leftmost/rightmost point of the corresponding circle, it is 
again split at that point. 

Step-2 Sweep a vertical line from left to right to compute the cells of the arrangement A{C). We also 
compute the ID field of each cell during the sweep. 

Step-1 

Consider a circle Ci € C. Each circle Cj £ C, j ^ i, creates an arc along the boundary of Ci that indicates 
the portion of the boundary of Cj that is inside Cj. In order to split Cj into pseudo segments, we need to 
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compute the minimum mimbcr of rays from the center of Ci that are required to pierce all the arcs a*, 
j = 1,2, ... ,n, j ^ i. This can be computed using the 0{n) time algorithm for computing the minimum 
geometric clique cover of the circular arc graph provided the end-points of the circular arcs are sorted [14]. 
But, we need to sort the end-points of the circular arcs along the boundary of Ci. Thus, the splitting of 
all the circles in C into pseudo segments need O(n^logn) time. Tamaki and Tokuyama [26] showed that 
the number of pseudo segments may be 0{ni) in the worst case. Recently Aronov and Sharir [3] showed 
that number of pseudo-segments generated from n unequal circles is at most n2+^, where e can be made 
arbitrarily small. 

Step-2 

In this step, we sweep a vertical line from left to right exactly as in [17] to compute the arrangement A{C). 

During the sweep, four types of events may occur: (i) leftmost point of a circle (ii) rightmost point of a circle, 

(iii) an end-point of a pseudo-segment that is not of type (i) or type (ii), and (iv) intersection point of two 
pseudo-segments. The events of type (i), (ii) and (iii) are initially inserted in a heap H. The events of type 

(iv) are inserted in H when these are observed during the sweep. The sweep line status data structure C 
stores the edges intersected by the sweep line at the current instant of time. Each pair of consecutive edges 
indicate a cell intersected by the sweep line. Each cell rj intersected by the sweep line (indicated by a pair 
of consecutive edges in the sweep line status) is attached with a balanced binary search tree containing 
the radii of the circles overlapping on that cell. Each time an event having minimum x-coordinate is chosen 
from H for the processing. The actions taken for each type of event is listed below. 

• While processing a type (i) event corresponding to a circle C, a new cell rj and two new edges, say ei and 

62, of A{C) take birth. These two new consecutive edges are inserted in C. If the new cell rj arrives inside 
an existing cell 7]' in the sweep line status £, then r^, attached to the cell 77, is created by copying Trj> 
and inserting the radius of the circle C in it. The ID field attached with the cell r] is the largest element 
in Tr,. 

• Type (ii) events are also handeled in a similar fashion. Here two edges are deleted from the sweep line 

status jC. Thus, a cell will also disappear from £. 

• At a type (iii) event one edge leaves the sweep-line and a new edge appears on the sweep line. Here, 

excepting this change on the sweep line, no other action is needed. 

• While processing a type (iv) event, an old cell rj' disappears from the sweep line and a new cell rj takes 

birth. If the event is generated due to the intersection of edges ei and 62 corresponding to the circles C'l 
and C2, then r,, is obtained by doing an O(logn) time updating of r^/. If r] is inside (resp. outside) of 
the circle Ci, then the radius of Ci is inserted in (resp. deleted from) r^/ to get r^. Finally, the largest 
element of r^, is attached as the ID field of the cell 77. 

Since the number of type (i) events is 0(n), and each type (i) event needs 0{n) time (for copying a heap for 
the new cell), the time needed for processing all type (i) events is O(n^) in the worst case. The number of 
type (ii) and type (iii) events are 0{n) and 0{ni) respectively. As mentioned above, processing each type 
(ii)/type (ii) event needs 0(1) time. The number of type (iv) events is 0{nP) in the worst case, and their 
processing needs 0{n^ logn) time. The point location in the arrangement A{C) of pseudo-segments is similar 
to that in the arrangement of line-segments. Using trapezoidal decomposition of cells, one can perform the 
query in O(logn) time [23]. Thus, we have the following result: 

Theorem 2. Given a set of circles of arbitrary radii, the preprocessing time and space complexity of the LCQ 
data structure are O(n^logn) and O(n^) respectively, and given an arbitrary query point, the largest circle 
containing it can be reported in O(logn) time. 
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3 QMEC problem for convex polygon 



Let P be a convex polygon and {pi,P2, ■ ■ ■ ,Pn} be its vertices in anticlockwise order. The objective is to 
preprocess P such that given an arbitrary query point q, the largest circle Cq that contains q but not 
intersected by the boundary of P can be reported efficiently. Needless to say, if q lies outside or on the 
boundary of P, Cq is a circle of infinite radius passing through q. So, the interesting problem is the case 
where q lies inside P. Needless to mention that here Cq is an MEC inside P. The medial axis M of P is 
the locus of the centers of all the MECs' inside P. Let c be the center of the largest MEC inside P^ (see 
Figure 1(a)). The medial axis consists of straight line segments and can be viewed as a tree rooted at c [10]. 
To avoid the confusion with the vertices of the polygon, we call the vertices of M as nodes. Note that, the 
leaf-nodes of M are the vertices of P. Let us denote an MEC of P centered at a point x G M as MEC^ and 
let Ax be the area of MECx- 




Observation 1 As the point x moves from c along the medial axis towards any vertex pi € P (leaf node of 
M), Ax decreases monotonically (see Figure 1(a)). 

Proof. Follows from the convexity of the polygon P. □ 

The medial axis M partitions the polygon P into n convex sub-polygons such that each sub-polygon Pi 
consists of a polygonal edge PiPi+i and two convex chains of M , one starting at pi and other starting at Pi+i 
(see Figure 1(b)). This partitioning can be achieved in 0(n) time since M can be computed in linear time 
[10]. Moreover, M can be preprocessed in 0{n) time so that the sub-polygon containing any query point q 
can be located in O(logn) time [18]. 

Lemma 1. The polygon P can be partitioned in 0{n) time such that given any arbitrary query point q, the 
edge of M closest to q can be reported in O(logn) time. 

Proof. We consider each Pj separately, and compute the medial axis of the convex chain from pi to pi+i (a 
portion of M). This needs 0(/Lij) time [1], where /itj is the number of nodes in M that appear as the vertices 
of Pi. Thus, for the entire polygon P, the total time complexity is 0(X]i=i Mj) = 0{n), since the number of 
edges of M is 0(n), and each edge of M appears in exactly two sub-polygons. If the query point q appears 
in Pj, then we can locate the edge of M that is closest to q in 0(log/Ui) time using point location in planar 
subdivision [18]. □ 

® There can be infinitely many such MECs of equal radius inside P in a degenerate case. But for the sake of simplicity 
we will assume that the largest MEC inside P is unique 
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(a) (b) 
Fig. 2. (a) Partitioning the medial axis M, and (b) The subtree M"^ for a query point q 



Now wc will describe how to solve the QMEC problem for a convex polygon. Assume that the query point q 
lies inside the sub-polygon Pi, that is incident to the edge PiPi+i of P. Let c' denote the center of the largest 
MEC containing q. Note that, c' will lie either on the path from pi to c (denoted by pi ~ c) or on the path 
from to c {pi+i ~ c) on M. Let us assume that c' lie on the path pi ^ c. We use Lemma 1 to identify a 
point X on the path pi ^ c that is closest to q in O(logn) time. The MECx must contain q. 

By Observation 1, we can locate c' by performing a binary search on the path c ~ x that finds two consecutive 
nodes v and v' on the path such that MECy encloses q, but MECyi does not. In degenerate case v may be 
X and v' is its previous node on the path c ^ x. Since the path lies on a tree representing the medial axis 
M, we can use level-ancestor queries [6] for this purpose. After computing v and v' , the exact location of c' 
can be determined in 0(1) time. Thus, we have the following theorem: 

Theorem 3. A convex polygon on n-vertices can be preprocessed in 0{n) time and space so that the QMEC 
queries can be answered in O(logn) time. 

4 QMEC problem for simple polygon 

Our approach for solving the QMEC problem in a simple polygon P is based on the divide and conquer 
strategy, and it uses the tree structure of the medial axis M. Here again the leaf nodes correspond to the 
vertices of the polygon. The internal nodes correspond to the points on M such that the MEC centered at 
eac;h of those points touches 3 or more distinct points on the boundary of P; We use J\f to denote the set of 

internal nodes of M . 

For the sake of simplicity in analyzing the algorithm, we assume that the MECs centered at the internal 
nodes of M have distinct radii. A point x G M, that is not a leaf, is said to be a valley point if for a sufficiently 
small S > 0, the MECs centered at points in M within a distance 5 from x are at least as large as MECx- 
We can similarly define the peaks in M. We assume that the number of peaks and valley points are finite. 
We use and to denote the set of valleys and peaks respectively. It is easy to observe that $ (lAf = ^, 
but C TV. 

Finally, we define a mountain to be a maximal subtree of M that does not contain any valley point except 
its leaves. Notice that, 

(i) Each mountain has exactly one peak. 

(ii) Each valley point is common to exactly two mountains, and it is a leaf for both the mountains. 

(iii) If a point x proceeds from a valley point of a mountain toward its peak, the size of MEC^ increases. 

Thus, if we partition M by cutting the tree at all the valley points, we get a set of mountains M = 
{Mi,M2,...,M\M\} (See Figure 2(a)) . 

We also need to consider another way of splitting the tree M as stated in Lemma 2. This aids in designing 
a data structure T for the query algorithm. 
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Lemma 2. [16] Every tree T with n nodes has at least one node n whose removal splits the tree into subtrees 
with at most \^~\ nodes. The node n is called the centroid ofT. 

Lemma 3. // q is the query point and M'^ is the maximal portion of M such that MECs centered in any 
point on enclose the query point q, then M'^ is a connected subtree of M (see Figure 2(b)). 

Lemma 4. // q falls outside the MEC^^, then is contained entirely in one of the subtrees obtained by 

deleting tt from M . 

Proof. Follows from the connectedness of (see Lemma 3). □ 

Lemmata 3 and 4 lead to the following divide and conquer algorithm for the QMEC problem for a simple 

polygon. 

In the preprocessing phase, we first compute the medial axis M. Next, we create a tree T whose root node 
is the centroid tt of M. The children of the root node in T are the centroids of the subtrees obtained by 

deleting tt from M. The process continues up to the leaf level. Each node f e T is attached with MECy. 
Note that, the MEC attached to the root node of T may not be the largest MEC in P. 

During the query with a point q, we need to consider two cases: (i) q lies inside the MECy for some vertex 

V of the medial axis M, and (ii) q docs not lie in the MEC of any vertex of the medial axis. We describe the 
method of computing Cq in Case (i). In Case (ii) (i.e., where Case (i) fails), then we identify the mountain 
Mi in which q lies. Next, we find Cq in Mi using the same method as in the Convex polygon Case, described 
in Section 3. 

The method of solving Case (i) is as follows. We test whether q lies in the MEC attached to the root node of 
T. If so, we report the largest MEC Cq using a data structure Query-in- Circle (or QiC in short), described 
below. If q does not lie in the MEC corresponding to the root node, then by Lemma 4, we need to search 
one of the subtrees of the root node. The search process continues until a node v' of T is identified such that 

MECy' contains q. 

During the search in the tree T, suppose we have identified a node v such that q lies inside MECy. Thus, 

V lies on the subtree M'. Here two important things need to be noted: (i) Cq may not be equal to MECy-, 
it may be some other MEC of larger area centered on M*, and (ii) M'' may consist of several mountains. 
The task of the QiC data structure attached to a node v of T is to identify the appropriate mountains in 

for searching the center of Cq. We also need another data structure, called MEC-in-Mountain (or MiM 
in short) that can report the largest MEC containing q with center on a given mountain Adi € M., provided 
MjflM^ 7^ 0. We now explain MiM and QiC procedures in detail, and then the divide and conquer procedure. 

4.1 MiM query 

Here we are given the polygon P and a mountain M^; we need to report the largest MEC centered at a point 
on Mi provided M, n ^ 0. Note that, if the center moves from any point x e fi to the peak of 
Mi, the MECs' arc strictly increasing. Thus, we can apply the algorithm proposed in Section 3 to identify 
the largest MEC containing g, and centered on Mj n M'^. The preprocessing time and space complexities 
are both 0(|Mi|), and the query time is 0(log |Mj|), where \Mi\ denotes the number of sides of the simple 
polygon that induces the edges of Mj. 

4.2 QiC query 

Here we want to solve a subproblem in which we know that the query point q falls inside an MEC centered 
at a given point v € M. We arc to prcprocess this information. In the query phase, given a query point 
q € MECy, we are required to report Cq, the largest MEC containing q. This problem is quite challenging 
since the locus M* of the center of Cq for possible choices of q satisfying above, is a subtree of M, and it 
may span several mountains. 
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Medial Axis path from p to c 

(a) (b) 
Fig. 3. (a) Bounding |<S|, and (b) Illustration on the number of MECs in .S^ that enclose a point x 



During the breadth-first search in 7", suppose we have already identified a vertex f in 7" such that MECy 
contains the query point q. But Cq may be some other MEC of larger area. We need to identify Cq. By 
Lemma 3, both the center of Cq and the node w of T are guaranteed to be on M^. 

Let TZv be the set of radii of MECs centered at the internal nodes of the subtree M* rooted at v, sorted in 
increasing order. 

Definition 1. An MEC C is called a guiding MEC corresponding to a node v of M* if 

• its radius is in TZy, 

• every MEC in the path from v to the center of C (both inclusive) is no larger than C, and 

• C overlaps with MECy. 

Let S be the set of all guiding MECs of the node v. Note that, a member in <S may be centered at the nodes 
as well as on the edges on M*. 

Preprocessing steps in QiC We perform the following steps in the preprocessing phase to compute <S 
attached to a node v oi M. Let M* be the subtree of M attached to node v. 

1. Perform a breadth first search in M* starting at v, and it recursively proceeds as follows: 

At each step (at a node v' G M*), if MECy' does not overlap with MECy, the recursion stops along 
that path; otherwise, two distinct cases need to be considered. We check whether the radii of MECy' 
and MECy" are consecutive elements in TZy, where v" is the predecessor node of v' in M*. 

• If so, put MECy' in S and recursively explore all the paths incident at v'. 

• Otherwise, compute all the MECs with center on the line segment {v, v') whose radius matches with 
the elements in the array TZy, put them in <S, insert those points on {v, v') as the (dummy) nodes in 
the tree M, and then recursively explore all the paths incident at v' in M*. 

2. Attach the mountain-id with each C G S. This is available while performing the breath-first search. This 
will allow us to invoke the MiM query for a particular mountain. 

3. Attach each circle in <S with the corresponding mountain in M. 

4. Create a LCQ data structure with the circles in <S, and attach it with node v. 

Lemma 5. For any r G TZy, the number of circles in S of radius r attached with node v is bounded by a 
constant. Furthermore, S can be computed in 0(|i?„|) time. 
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Proof. Consider any r e TZy Let S,. he the; MECs of radius r in S. It sufSccs to show that \Sr\ is bounded 
by a constant. For convenience, let us assume that Sr does not contain a MEC centered at a node of M. 
This will not affect us because we have assumed that MECs centered at nodes of M have distinct radii, so 
at most one MEC in Sr can be centered at a node. 

Let r„ be the radius of MECy of the root node v of M*. Clearly, < r (by Definition 1). Also recall that 
every MEC in <S must (at least tangentially) intersect MECy. See Figure 3(a) for an illustration. Therefore, 
every MEC in Sr must lie entirely within a circle x of radius r„ + 2r centered at v. Thus, we need to prove 
that the number of guiding circles of radius r at node v inside x is bounded by a constant. 

Let us consider a point x & P . Let <S^ C ,5^ be a set of MECs that enclose x. Let C be any MEC in <S^f and 
c be its center. Let pi and p2 be the two points at which C touches the boundary of the polygon P. The 
chord [pi,P2] must intersect the medial axis (see Figure 3(b)). Note that, the points p and c lie in the two 
different sides of [pi,P2]- On the contrary, if p and c lie in the same side of [p'i,P2], where p'l and p'2 are the 
points of contact of the said MEC and the polygon P, then we can increase the size of the MEC by moving 
its center c towards p along the medial axis (see Figure 3(b)). Thus, C ^ Sr. Thus, we have Zpixp2 > tt/2. 
These angles subtended by the MECs in <S^ are disjoint implying that \Sr\ < 4. In other words, any point 
inside the circle x can be enclosed by at most four different circles from Sr- We need to compute \Sr\. Let 
us consider a function f{x) = number of circles in Sy that overlaps at the point x, x € x- fi^) ^ 4 for all 
X G X. The total number of circles in Sy can be obtained as follows: 

Total area of circles in Sr < /^^ ^j^^ f{x) dx dy < 47r(r„ + 2r)^. 

Therefore, |5J < Mi^^+M! < MM! = gg. 
Thus, the first part of the lemma is proved. 

The time complexity follows from the fact that the breadth first search in M* needs 0(|M*|) time. The time 
for computing the members in S is ^ren„ \^r\ = 0{\TZy\) (by the first part of this lemma). □ 

Query algorithm in QiC Given a query point q, we first traverse the tree T to identify a node v such that 

q e MECy. Note that, the MECs at the nodes of the subtree rooted at v may contain q; but the MECs 
corresponding to all other nodes in T will not contain q (see Lemma 4) . Let S be the guiding circles attached 
with node v, p GTZy be the radius of the largest guiding circle in S that contains q, and Sg be a subset of <S 
that has radius p and contains q. Sq can be obtained from the LCQ data structure attached with node v. By 
Lemma 5, we have l^gl < 36. In order to report Cq, we need the following: 

Step 1: an algorithm to identify the mountain associated with each circle in Sq, 

Step 2: to locate the largest MEC containing q in each of these mountains using the MiM query algorithm, 
Step 3: to report the largest one among the MECs' obtained in Step 2 as Cq. 

We first devise an algorithm for Step 1. The necessary algorithm for Step 2 is already available in Subsection 
4.1. We then prove the necessary result to ensure the statement stated in Step 3. 

Algorithm for Step 1: 

Consider a path 11 from w to a leaf of M*, and observe the size of the MECs'. Figure 4 demonstrates a curve 
f{t) where t denotes the distance of a point from v on the path 11, and f{t) denotes the radius of the MEC 
centered at that point. The guiding circles along the path U correspond to a subsequence of vertices along 
that path whose corresponding MECs' are increasing in size. 

Lemma 6. The guiding circles along a path U from v to a leaf of M* containing the query point q appear 

consecutively along 7J. 

Proof. Follows from the connectedness of (see Lemma 3). □ 
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Fig. 4. Proof of Lemma 7 

Consider tlie MECs' attached to the nodes in 77. Let v' be such a node whose corresponding MEC is largest 
among those containing q. Let M' (g M) be the mountain in which v' lies. Here two cases need to be 
considered: (i) v' is the peak of M' , and (ii) v' is not the peak of M'. In Case (i), we have already got the 
largest MEC centered on the path 77 and containing q. In Case (ii), we need to invoke MiM query algorithm 
to find the largest MEC centered on the mountain M'. 

Correctness of QiC 

Lemma 7. At least one of the circles in Sq is centered in the mountain in which Cq is centered. 

Proof. Since is a continuous subtree of M , if we explore all the paths in M* from node v towards its 
leaves, Cq is reached in one of such paths, say 77, and v' be a node on 77 such that the guiding circle MECyi 
is largest among those which contain q. Note that, any point on the path 77 closer to v than v' can not be 
the center of a larger MEC (see Definition 1). Let the center Cq of Cq be a point on 77 that is in a different 
mountain to the right of v' . Here again two situations need to be considered: (i) the function f{t) increases 
monotonically from v' to c,, and (ii) the function f{t) from v' to Cq is not monotonic. In Case (i) v' and Cq 
lie in the same mountain. In Case (ii), between v' and Cq there is a point a on the path 77, such that the 
radius of MECa is less than that of MECyi . Also, there exists another point /? on the path 77 between a 
and Cq such that the radius of MECp is equal to that of MECyi . Since the radius of MECp matches with 
an element of TZ, MECp is also a guiding circle. Moreover, from the continuity of M', the MEC centered 
at (3 must contain q. So, if Cq does not lie in the mountain of v' , it must lie in the mountain containing /3. 
Thus, the lemma follows. □ 

Lemma 8. The preprocessing time and space complexities for the QiC query are 0(|M|log^|M|) andO{\M\log\M\) 
respectively. Queries can be answered in 0(log^ |M|) time. 

Proof. Compuing S requires 0{\M* \ log |M* |) time because we need to sort the elements in 7?^. The members 
in S can be stored in LCQ data structure in 0(|7?.„| log^ \TZv\) time and 0(|7?.„| log \TZv\) space (see Theorem 
1) and queries can be answered in log^ l-^^ |- 

In the query phase with a query point q, we identify a constant number of guiding circles Sq attached to 

node V that contains q. Next, we call MiM queries in their associated mountains; this takes log \M* \ time (see 
Theorem 3). By Lemma 7, the result of one of the MiM queries will be the largest MEC with center on M* 
that contains q. Thus the query time complexity follows. □ 

4.3 Query algorithm for finding Cq 

Before we start the divide and conquer, we compute the set Ai of mountains and preprocess each of them 
for MiM query. Since is a partition of the medial axis M, all the mountains can be preprocessed for MiM 
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query in 0(|M|) time. In Lemma 8, it is shown that the total preprocessing time needed for the QiC queries 
at every node of T is 0{\M\ log^ \M\) using 0{\M\ log \M\) space. 

In the query phase, we start at the root level of T and check if q falls inside the MEC centered at the root. 
If yes, we find Cq using the query algorithm for QiC. This takes 0(log^ n) time (see Lemma 8). Otherwise, 
we proceed in the appropriate subtree of the root whose corresponding sub-polygon contains q. In order to 
choose this sub-polygon, we need another data structure as stated below. 




Fig. 5. The divide and conquer search structure 

Recall that, the medial axis M partitions P into n cells. Let tt be the centroid node that corresponds to the 
root of T. Let tt have k children. In other words, if we consider the AIECtt, it touches P at fc different points. 
This gives birth to k sub-polygons (as illustrated in Figure 5 with A: = 5). The centroid of each sub-polygon 
is a child of tt. We attach a first-level-tag i with each cell in the i-th sub-polygon, for i = 1,2, . . . , k. Next, 
we consider the children of vr (the nodes in the second level of T) in a breadth first manner. For each sub- 
polygon, consider its centroid. The MEC of that node again partition that sub-polygon into further parts. 
We attach a second-level-tag to each cell of that sub-polygon as we did for the root. After considering all the 
children of tt, we go to the third level, and do the same for attaching the third-level-tag to the partitions of 
P. Since the number of levels of T is O(logn) in the worst case, a cell of P may get O(logrt) tags. Thus each 
cell is attached with an array TAG of size 0{\ogn) containing tags of O(logn) levels. This needs 0(7ilogn) 
time and space in the worst case. 

By point location in the planar subdivision of P, we know in which partition Q of P the query point q 
lies. While searching in the tree T, if q does not lie in MECu of a node u in the i-th level, we choose 
the appropriate subtree of u by observing the i-th entry of the array TAG attached to the partition Q, and 
proceed in that direction. Thus, the overall query time complexity includes (i) O(logn) for the point location 
in the subdivision of P, (ii) O(logn) time for traversal in T, (iii) 0(log^n -I- K\ogn) time for identifying 
the largest guiding circles attached to node v (in its LCQ data structure) if node v is observed first during 
traversal of T, such that MECy contains q, and the MiM queries for K mountains if K circles are output of 
step (iii). In Lemma 5, it is proved that K is bounded by a constant. Thus, we have the following theorem. 

Theorem 4. A simple polygon can be preprocessed in 0(nlog'^ n) time using (nlog^ n) space and the QMEG 
queries can be answered in 0(log^ n) time. 



5 QMEC for Point Set 

The input consists of a set of points P — {pi,P2, . . . in M^. The objective is to preprocess P such that 
given any arbitrary query point q, the largest circle Cq that docs not contain any point in P but contains q, 
can be reported efficiently. Observe that, if q lies outside or on the boundary of the convex hull of P, we can 
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draw a circle of infinite radius passing through q. So, we shall consider the case where q lies in the proper 

interior of the convex hull of P. 

An MEC centered at a Voronoi vertex touches at least three points from P. We assume that the MECs 
centered at Voronoi vertices are of distinct sizes. For our purpose, we also compute some artificial vertices, 
one on each Voronoi edge that is a half line. We must compute these artificial vertices carefully to ensure 
that the following conditions hold. 

1. Every MEC centered at an artificial vertex must be larger than MECs centered at Voronoi vertices, and 

2. the MECs centered at artificial vertices should not overlap pairwise within the convex hull of P. Surely, 
they overlap outside the convex hull of P. 

The second condition ensures that there exists no query point q which can be enclosed by more than one 
MEC centered at artificial vertices. This second condition makes the choice of artificial vertices somewhat 
tricky, but it is a simple exercise to see that we can choose the artificial vertices in O(n^) time. We use the 
unqualified term vertex to refer cither to a Voronoi vertex or an artificial vertex. 

Now, consider the planar graph with both the Voronoi vertices and the artificial vertices. Let vhe & Voronoi 
vertex and let MECy be the MEC centered at v. A path Vy = {v^ = v,v'^ , . . . ,v'') from v in the graph is 
said to be a rising path with respect to v if 

• MECs centered at vertices other than and v'^ are strictly smaller than MECy, and 

• MECyk is strictly larger than MECy. Note that v'^ may be an artificial vertex, but the other vertices 
in the path Vy are surely Voronoi vertices. 

The last edge {v'^~^,v'^) is called a rising edge with respect to the vertex v. Since MECyk-i is smaller than 
MECy, but MECyk is larger than MECy, there is exactly one MEC centered on the edge {v''~^,v'') that 
equals the size of MECy. Let us denote this MEC by MEC-p^. If MEC-p^ overlaps MECy, then the rising 
edge; (i,'*'^^, w*^) is called an overlapping edge with respect to vertex v. Let Oy be the set of overlapping edges 
with respect to v. Observe that Oy, for a given vertex v, can be computed in 0(n) time via a breadth first 
search from v. The preprocessing and query procedures are given in Procedures 1 and 2. 



Procedure 1 Preprocessing steps 
1: INPUT: Set of points P in M^. 

2: Compute the MECs centered at (both Voronoi and artificial) vertices, and store them in LCQ data 
structure. 

3: For each vertex v, compute Oy using a breadth first search. 



We are now left with showing that Procedures 1 and 2 are correct and bound their complexities. We address 
the latter first. We begin with a lemma, which can be proved essentially using the proof of Lemma 5. 

Lemma 9. For any internal vertex v in the Voronoi diagram of P, \Oy\ is hounded by a constant. 

Proof. Let MECy be the circle centered on v. Consider any overlapping edge e = {vi,V2) in Oy. Assume 
without loss of generality that the MEC at V2 is larger than the MEC at vi . By definition, there is a point v' 
on e such that MECy' has the same radius as MECy and that MECy' intersects MECy. Let pi and P2 be 
the two points in P that touch the MEC at v'. The chord piP2 intersects the edge e somewhere between vi 
and v' (see Figure 3). Otherwise, the MEC at will not be the first MEC from vi to V2 that equals MECy 
in size. Therefore, we can use the same idea from Lemma 5 to bound the number of overlapping edges. □ 
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Procedure 2 Query steps 

1: INPUT: a query point q along with the LCQ data structure containing MECs centered at vertices, and 

Oy for every internal vertex v in the the Voronoi diagram of P. 

2: Find the largest MEC C in the LCQ data structure containing q. Let v'^ be the center of C. 

3: if v'^ is an artificial vertex then 

4: Report Cq = C as the largest circle containing q that is centered on the edge containing w^. 

5: Exit 

6: end if 

7: C^C 

8: for all edges e e Oyq do 

9: if MECs centered on e do not enclose q then 

10: Continue to next edge in O^q 

11: end if 

12: Let be the largest circle centered on e that encloses q. 

13: C<«-max(C,C"=) 

14: end for 

15: Report C 

Lemma 10. Given that we can construct the LCQ data structure forn circles in 0{p{n)) preprocessing time 
with a space complexity of 0{s{n)) and queries answered in 0{q{n)) time, Procedure 1 (for preprocessing) 
takes 0{p{n) + n^) time and 0{s{n) + n) space, and Procedure 2 (for query answering) takes 0{q{n)) time. 

Proof. The complexity bounds for LCQ data structiirc arc added for the obvious reason that we use the LCQ 
data structure for creating and storing 0{n) MECs centered at the internal vertices of the Voronoi diagram. 
Line number 3 of Procedure 1 performs 0{n) breadth first searches, hence we added an O(n^) term to the 
preprocessing time. As a consequence of Lemma 9, our space requirements is limited to 0{n) and, more 
importantly, the query time does not incur anything more than q{n). □ 

Lemma 11. Consider any cycle H in the Voronoi diagram of P. Let C'h be any MEC centered at some 
point on H. Then, there exists another MEC C'^ centered at some other point on H that does not properly 
overlap Ch- 

Proof. Clearly, any cycle in the Voronoi diagram of P must contain at least one point from P inside it. Let 
p e P be such a point that lies inside the cycle H (see Figure 6). Let Ch be any MEC centered at some 
point on H] ch be its center. Consider the line connecting ch and p. It intersects H at another point c'^. It 
is easy to see that the MEC C^, centered at c'^, will not properly overlap with Ch. Because, in that case p 
will be properly contained within Ch and C^. □ 

Lemma 12 (Unique Path Lemma). Let C and C be any two distinct but overlapping MECs with center 
at c and respectively. There is exactly one path from c to c' along the Voronoi edges such that every MEC 
centered on that path encloses C (iC . 

Proof. The structure of the proof is as follows. We provide a procedure that constructs a path iT(c, c.) from 

c to c' along the Voronoi edges, and ensure that every MEC centered on that path encloses C O C . As a 
consequence of Lemma 11, the path does not form an intermediate cycle and terminates at c'. Finally, we 
again use Lemma 11 to show that no path V, other than n{c,c'), exists between c and c' such that every 
MEC centered on V contains C OC . Throughout this proof, we closely follow Figure 7 in order to keep the 
arguments intuitive. To keep arguments simple, we assume that c and c' are Voronoi vertices. When c and 
d are not Voronoi vertices, then also the same argument follows. 

Let a be the number of points in P that C touches. These a points partition C into a arcs. The degree 
of the corresponding Voronoi vertex c is also a because each adjacent pair of points from P that lie on 
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Fig. 6. Illustration of Lemma 11; the edges in H are shown darker 

the boundary of C will induce a Voronoi edge incident on c and vice versa. These Voronoi edges and their 
corresponding arcs are denoted by e;^ and sPq, for 1 < j < a. 

Consider the other MEC C C and centered at a vertex c') that overlaps with C. C intersects C at two 
points ti and t2] both ti and t2 must lie in one of the a arcs of C (due to the emptiness of C"). Let us name 
this arc by s;^. Consider the edge e'^ = (c, C2) that corresponds to the arc s;^. The other end of e^, i.e., the 
vertex C2, is called the next step from c toward c' and denote it as ns(c, c'). Consider the following code that 
generates a path denoted by 11 (c, c'): 



Procedure 3 -/T(c, c') Computation 

1: n(c,c') <(- (c) 
2: next ^ c 
3: repeat 

4: next ns( next, (/) 
5: Append next to n(c, c') 

6: until next equals c' {This is the only terminating condition.} 



Let n(c, c') as constructed above be (ci = c, C2, . . . , Ci, Q+i, . . . , c'). Let C2 denote the MEC centered at C2. 
If C2 is the circle C", then the procedure terminates and, as required, every MEC in the edge (c, C2) encloses 

CnCs = CnC". 

Therefore, consider the case where C2 is not C . Let pi,p2 € P be the points at which C and C2 intersect; 
Pi,P2 are the end points of the arc s;^ that defines the next step move toward d (in Figure 7, j is 2). 
Therefore, by definition, ti and i2 lie on the arc s^. Notice that CnC" (shown shaded in Figure 7) is shaped 
like a rugby ball with ti and ^2 at its end-points. One side of C n C" (called the initial side) is in C and the 
other side (called the final side) is in C . Clearly, ti and t2 are inside (or on the boundary of) every MEC 
centered on the edge e^. Otherwise, as we go from C to C2, there will be a circle that touches the final side 
of C n C", but that would mean that we have either 

— reached C", which contradicts our assumption that C2 is not C", 

— or found a MEC that contains C , which contradicts the fact that C is itself an MEC. 
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Fig. 7. Illustration of n{c,c'). 



We now make two observations: (i) C touches the initial side, but (ii) no other MECs centered on (and C2 
in particular) touches the final side Suppose (for the sake of contradiction) that there is a MEC C* centered 
on Cq that touches the final side of C fl C" at, say, some point t*. It is easy to see that C* will contain C 

becaiise it touches f* on C" and contains ti and t2, which are also on C" this is a contradiction due to the 
observations (i) and (ii), stated above. Thus, it is clear that C fl C" is properly contained within Cfl C2. 

Consider two adjacent vertices Cj and Ci+i along iT(c, c') with MECs Ci and Ci+i centered on them, respec- 
tively. The above argument can be easily extended to give us the following: 

Ci n C c Ci+i n C. 

Therefore, we can conclude that every MEC along n{c,c') encloses C fl C . Given Lemma 11, we can also 
conclude that n{c,d) does not form a cycle. The only stopping condition is when we actually reach c', so 
iT(c, c') terminates at c' in 0{n) steps. Therefore, -/7(c, c') fulfils our requirements. 

To complete the proof of this lemma, we must show that -/7(c, c') is the only required path. For the sake of 
contradiction, assume that there is another path 77' such that every MEC centered on 77' contains C n C". 
Then, there are two distinct paths from c to c' such that every MEC centered on both paths overlapped with 
C n C". Clearly, there must be a cycle when the two paths are combined. From Lemma 11, we know that 
there are pairs of MECs in the cycle that will not overlap each other. This is a contradiction, thus proving 
that 77(ci, C2) is the only required path and concluding the proof of the lemma. □ 

To illustrate the usefulness of the unique path lemma, consider the example depicted in Figure 8. Let 
(centered at d^) be the MEC returned by the LCQ data structure when queried with q. Suppose that Cq 
(centered on Cq) is the largest empty circle containing q. Suppose Procedure 2 reports C'q centered on c'q, 
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Fig. 8. Illustrates the usefulness of the Unique Path Lemma. 

which lies on an overlapping edge (a, 6). In terms of size, let < C'g < Cq. Such a behavior will render our 
algorithm incorrect. However, this incorrect behavior is only possible when Cq does not lie on an overlapping 
edge with respect to c^, and therefore. Procedure 2 may never find it. This may happen in two possible 
situations: 

Situation 1: The path shown in thick continuous segments is 77 (c',Cq). Since Procedure 2 reported 
c'q, the edge (a, 6) is an overlapping edge, the MEC at b must be larger than C*. However, the LCQ 
data structure did not report the MEC at b, but rather reported C. While q gC^ and q G Cq, q is not 
contained within the MEC at b. Prom the unique path lemma, clearly, b cannot be in n{c'^,Cq). Thus 
such a situation is impossible. 

Situation 2: Some other path that does not contain (shown using dashed line segments) is 
n{cfl,Cq). This case is possible and the edge (c, &), shown in Figure 8, is an overlapping edge through 
that dashed path. Therefore, the edge (c, b) will be considered by Procedure 2 and Cq will be reported 
correctly. 

Lemma 13. If a set P of points from are preproeessed by Procedure 1, then Procedure 2, when invoked 
with a query point q, correctly reports the largest circle that contains q but is devoid of points from P. 

Proof. Let Cq be the largest circle that contains q but is devoid of points from P. In short. Procedure 2 must 
report Cq. If Cq is centered on a Voronoi vertex, then, clearly. Procedure 2 reports it correctly. 

Wc now show that if q is contained by a MEC centered at an artificial vertex v, then also the algorithm 
reports the correct Cq. Let e be the half line edge containing v. Procedure 2 only searches edge e, which, we 
claim is suflacient. Suppose for the sake of contradiction that Cq is centered on some other edge e'. Clearly, 
e' cannot be a bounded Voronoi edge, because the MEC at v is larger than MECs centered on bounded 
Voronoi edges. Recall that our construction of artificial vertices ensures that no two MECs centered on 
artificial vertices will overlap inside the convex hull of P. Therefore, e' cannot be an edge that is a half line 
either, because, the MEC at v contains q, so the MEC centered at the artificial vertex on e' cannot contain 
q. Therefore, such an e' cannot exist and we can conclude that the algorithm correctly reports the largest 
MEC centered on some point in e that contains q as Cq. 

For the rest of the proof, we assume that Cq is not centered on a Voronoi vertex and q is not enclosed by any 
MEC centered at an artificial vertex. Recall from Procedure 2 that C is the largest MEC that is centered on 
a vertex and contains q. Clearly, C H 7^ as it at least contains q. From Lemma 12, there is exactly one 
path n{c'^, Cq) from the center of to the center of Cq such that every circle in that path contains C n Cq. 
Clearly, MECs centered at vertices in n{c'^,Cq) other than the centers of C and Cq are smaller than C; 
otherwise, the LCQ data structure would not have chosen C^. However, Cq is larger than C. Consider the 
two vertices connected by the edge that contains the center of Cq. The MEC centered on one of them must 
be strictly larger than C, while the MEC on the other must be strictly smaller than C^. Therefore, it is 
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easy to see that Cq is centered on an overlapping edge. Since the algorithm searches through all overlapping 
edges, it will find and report Cq correctly. □ 



Lemma 10 coupled with the line sweep method of implementing the LCQ data structure and Lemma 13 
immediately lead to the following theorem. 

Theorem 5. Given a set P of points in M^, we can preprocess P in O(n^logn) time and 0{n'^) space so 
that the resulting data structure can be queried for the largest empty circle containing the query point q in 
O(logn) time. 

6 QMER problem 

The input consists of a set of points P = {pi,p2, . . . ,p„} in a rectangular region A. An axis-parallel rectangle 
inside A is said to be an empty rectangle if it does not contain any point of P. An empty rectangle is called 
maximal empty rectangle (MER) if no other empty rectangle in A properly contains it. Here our objective 
is to preprocess P such that given any arbitrary query point q € A, the largest area rectangle L(q) inside 
A that does not contain any point in P but contains q, can be reported efficiently. First observe that L{q) 
is an MER. Let M denote the set of all possible MERs in A. In [20], it is shown that \M\ = 0{n'^) in the 
worst case. 

In the preprocessing phase, we partition A into a set C of cells, such that for every point q inside a cell 
cG C, the largest MER containing q is the same. This is achieved by drawing horizontal and a vertical lines 
through each point in P. This splits A into the set C of O(n^) cells. Observe that for any cell c E C and 
any MER M„ € M, either c fl = c or c fl = (p. Therefore, if the query point q lies inside c G C, we 
need to report the largest MER containing c. We choose a representative point inside each cell c G C. Let 
Q be the set of representative points. We compute all the MERs using the algorithm in [20], and sort them 
with respect to their area. We also construct an augmented dynamic range tree T with the points in Q in 
0(n^ logn) time and space [19]. Next, we process the members in M in order. For each G M, we identify 
the set Qa of points in Q that are inside Ma ■ We store a pointer to Ma along with each point in Qa and 
then delete Qa from T. This step takes 0{L + logn) time [19], where L is the number of points inside M^. 
After processing all the MERs in M, we have stored a pointer to the largest MER along with each q G Q. 
Thus, we have the following theorem: 

Theorem 6. A set of n points in a rectangle A can be preprocessed in O(n^logn) time and space so that 
the largest empty rectangle query containing the query point can be answered in O(logn) time. 

Moreover, it is not hard to see that we can construct examples, where there are J?(n^) cells, so that the MER 
containing each of these cells is combinatorially different, i.e., the boundary of any of the two MERs are not 
incident to the same set of points in P. This suggests that in order to answer queries in polylogarithmic 
time, we need to somehow store J7(n^) cells in a data structure, and hence it is very unlikely to improve the 
cost of preprocessing in order to maintain O(logn) query time. 

7 Future Work 

Our focus in this paper has been in terms of understanding which problems can be solved within subquadratic 
preprocessing time, while maintaining the polylogarithmic query time. At this stage the central problem here 
is to understand whether the preprocessing time for the QMEC problem for the point set case can be tightened 
to a subquadratic bound. A possible lead for improvement is as follows. If we were to insert q into the set of 
points P, and compute the Dehmnay triangiilation of the new set, then the query circle, Cq, is a circumcircle 
of one of the triangles t incident to g - in fact the angle subtended at g in f will be greater than 7r/2. The 
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running time for inserting q in the Dclaunay triangiilation of P is proportional to the degrcie of q. One may 
be tempted to use a (randomized) incremental algorithm for constructing Delaunay triangulation, but there 
are cases in which the degree of q can be linear. This approach may lead to a practical and a simpler way to 
handle these types of queries, but this remains to be seen. Alternatively, one may look into the divide and 
conquer algorithms for computing the Voronoi diagram, and see whether in the "merge step" , the maximal 
empty circles can be maintained. Alternatively, one may try to use the planar separator theorem to partition 
the Voronoi diagram, recursively, and for the separator vertices (point), biiild an appropriate structure, so 
that preprocessing can be performed in subquadratic time, and the queries can be answered in sublinear 
time. 

It will be desirable to improve the preprocessing cost in the case of simple polygons. A real challenge will be 
to match the complexity in this case to exactly that of the convex polygon case. 

While we have studied a few canonical problems, there are several other variants that are as yet untouched. 
We can also ask similar questions on multidimensional geometric sets, but we suspect that the curse of 
dimensionality might restrict us to approximations. 
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